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1. Summary 


*An Euler flow solver has been developed for predicting the airframe/propulsion 
' - integration effects for an aft-mounted turboprop transport. This solver employs a 

highly efficient multigrid scheme, with a successive mesh-refinement procedure to 
accelerate the convergence of the solution. A new dissipation model has also been 
implemented to render solutions that are grid insensitive. The propeller power 
effects are simulated by the actuator disk concept. An embedded flow- solution 
method has been developed for predicting the detailed flow characteristics in the 
local vicinity of an aft-mounted propfan engine in the presence of a flowfield 
induced by a complete aircraft. Results from test case analysis are presented. A 
user’s guide for execution of computer programs, including format of various input 
files, sample job decks, and sample input files, is provided in an accompanying 
volume. 


2. Introduction 

There is a strong need for computational aerodynamics codes for the analysis of 
complete aircraft configurations. A method for airframe/engine integration must 
correctly include all the complexity of the flow. Computational methods can 
provide detailed knowledge of the flow and aid in the design of superior 
configurations in a cost-effective manner. In recent years, significant progress has 
been made in the development of numerical methods for general three-dimensional 
analysis. However, much research work remains to be done in the area of 
transonic, rotational flow analysis for complex airplane configurations with 
multiple lifting surfaces. The rotational flow analysis is necessary for analyzing a 
turboprop transport because of the swirl behind the propeller plane. .. 

- The objective of the present study is to develop a new computer program that 
can simulate the complex flowfield produced by an advanced aft-mounted propfan 
aircraft. This includes the development of an embedded flow solution method for 
predicting the detailed flow characteristics in the local vicinity of an aft-mounted 
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propfan engine in the presence of a flowfield induced by a complete aircraft. 


This program can be used to predict/evaluate the airframe/engine integration 
effects of a complete wing/body/tail/aft-mounted propfan configuration, as well as 
the effects of engine power on lateral stability (Ref. 1). This program can also be 
used as a technology basis for the development of an aerodynamic design 
capability for aft-mounted propfan aircraft by coupling the program with an 
automated iteration procedure developed at Langley (Ref. 2). 


3. Global Flowfield Solution 

A solution method that can predict the global flowfield about a complete aircraft 
configuration with aft-mounted engines has been developed. The global flowfield 
solution method includes two major steps. These are: a. the generation of the 
computational grid, and b. the development of Euler flow solver. The numerical 
validation of the global solution method is also discussed to show how the method 
works. 

a. Grid Generation 

The analysis of a complete airplane, such as that shown in Fig. 1, requires the 
use of a grid system that will yield accurate solutions at affordable cost using 
existing supercomputers. The unstructured-grid approach (Refs. 3, 4) offers great 
flexibility in treating complex geometry. However, the unstructured-grid flow 
solver requires a large central memory, and the solution algorithm is less efficient 
and less accurate than the structured-grid counterpart on a comparable grid. The 
rectangular-grid approach (Ref. 5) is an alternative, provided automatic mesh 
refinement can be done effectively and accurately. The present work utilizes a 
' structured surface-fitted grid. In this approach, the grid on the configuration 
surface is generated first, and then the flowfield grids are obtained using an elliptic 
grid generation scheme. This offer two major advantages: firstly, the structured- 
grid requires minimum storage for the flow solution algorithm since the grid points 
are well ordered, and secondly, the surface-fitted grid allows for easier and more 
accurate implementation of the boundary conditions compared to a non-surface- 
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fitted grid. 


To generate grids for a complex aircraft configuration, a block-structured grid is 
used. The complete flowfield is divided into a number of convenient blocks so that 
smooth 3D grids can be generated in each block. Fig. 2 illustrates the present 
blocking strategy. A surface passing through the middle of the nacelle together 
with a surface passing through the T-tail plane (shown in dashed lines) is used to 
divide the complete flowfield into four blocks. Within each block, a 3D elliptic 
grid generation technique is derived from the method of Thompson (Refs. 6, 7). 
Source terms based on boundary grid distributions are added to the grid generation 
equations to control grid spacing and grid orthogonality in the flowfield. Details 
on grid generation equations and the implementation of source terms are discussed 
in Appendix A. 

The grid topology used in the present work is a C-type grid about the wing, an 
H-type grid about the stmt, and the T-tail as illustrated in Fig. 3. The grid in the 
spanwise direction is of the H type. The crown and keel lines of the fuselage form 
one of the C lines. The leading edge and trailing edge of vertical tail are fitted with 
grid lines. The grid on the nacelle surface and plume that extends downstream to 
the far field is H type. To illustrate the grid topology used in the present analysis, 
the field grid on spanwise, streamwise and normal direction cuts are shown in Figs. 
4, 5 and 6 respectively. Fig. 7 shows a complete surface grid and part of the field 
grids for an aft-mounted propfan airplane. 

b. Euler Flow Solver 

The present Euler analysis method is based on a cell-centered, finite-volume, 
out-of-core multigrid Euler solver for a complete aircraft analysis (Ref. 1). The 
outermost iterative procedure carries out the successive mesh refinement. Starting 
' * •••=.. from the coarsest mesh level, the entire flowfield is divided into blocks and the 

entire flowfield data is stored on a solid state disk (SSD). The computation is done 
in a block-by-block manner through the use of highly efficient input/output data 
management. The flowfield in each block is advanced in time through one 
multigrid cycle (Ref. 8). Within each multigrid level, a five-stage, explicit, 
Runge-Kutta time integration scheme, together with an implicit residual smoothing 
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method, is used to update the flow variables to a new time level. The updated 
flowfield data are then moved to the SSD before another block of flowfield data is 
transferred to the central memory. The block boundary data, needed for the 
adjacent blocks in flux and dissipation term calculations, are saved in different 
location on SSD, so that they can easily be fetched during the calculation process.-^ 

At the present time, the program is written in such a way that the flowfield can 
be divided in an arbitrary number of blocks in both normal and spanwise directions 
to handle large problems. However, the general strategy is to keep the number of 
blocks to a minimum, such that the program can be executed efficiently in a vector 
machine. Optimum blocking strategy can be different to fit a specific machine 
architecture. For a machine with very large directly addressable memory (e. g., 
CRAY-2 with 128 MW of memory), the entire flowfield may be divided into one 
to three blocks depending on the geometry complexity of the configuration. The 
converged solution at a given mesh level for the entire flowfield is interpolated to 
the next finer mesh level to start the Euler calculation on that mesh. Successive 
mesh refinement stops when the solution has converged at the finest mesh level. 
The basic numerical scheme used is discussed in Appendix B. 

In the calculation process, the type of boundary condition along each boundary 
surface of a given block is specified. If it is an interblock boundary surface, 
conservation of mass, momentum, and energy are ensured with the adjacent block. 
Also, the flux term due to dissipation must be formulated precisely to satisfy the 
conservation laws across the block boundary. Due to the use of fourth-order 
dissipation terms, two layers of data from every adjacent block are needed in order 
to formulate a flux term that is identical to the one used in single-block 
calculations. For a far-field boundary surface, characteristic boundary conditions 
that take proper account of inflow and outflow conditions are implemented. Along 
a solid boundary surface, the conventional normal momentum equation is used to 
compute the surface pressure. The propeller power loading is simulated by an 
actuator disk where the total pressure, total temperature, and swirl distributions are 
prescribed (Ref. 1). 

Another way of simulating propeller power effects is to prescribed the thrust, 
normal force, and side force along the propeller disk and the work added to the 
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flow by the propeller. The flow variables downstream of the disk are related to 
their upstream values through the use of continuity, momentum, and energy 
equations. The major advantage of this method is that the effects of angle of 
attack, as well as the influence of side flow, can easily be simulated through the 
input of normal force and side force distributions along the propeller disk. A 
detailed formulation of this method can be found in Reference 9. Boundary 
condition methods for engine inlet and exhaust are given in Reference 10. For the 
global analysis of an aft-mounted propfan transport, some simplification has been 
made. The primary exhaust flow of the engine is treated as a solid body with the 
shape obtained from an isolated nacelle or propfan analysis. The flowfield inside 
the plume is not computed, but this is a small region compared with the propeller 
slipstream. Boundary condition implementations are reviewed in Appendix C. 

Important improvements have been made to the dissipation formulation to 
provide more accurate and reliable solutions for complex geometry. The spectral 
radius scaling of dissipation terms (Refs. 11, 12) has been adapted to the Euler 
code. The new formulation improves the consistency of the solution by decreasing 
sensitivity to the dissipation parameters and to the quality of the grid. This 
improvement over conventional dissipation formulations (Ref. 13) is especially 
pronounced on grids with a large aspect ratio or with nonsmooth mesh 
distributions (Ref. 14). 

The focus of this study is in the subsonic and transonic speed region. However, 
the analysis method developed under present study also works for supersonic flow 
and is directly applicable to High Speed Civil Transport (HSCT) type of 
configurations. Numerical validation of the method in the transonic region is given 
in the next subsection. 

c. Numerical Validation 

Employment of both successive mesh refinement and the multigrid acceleration 
scheme is effective in speeding up the convergence of the solutions to steady state. 
For a typical flow solution of a wing-body-winglet configuration containing 
200,000 grid points, three levels of mesh refinement are used in the calculations. 
' * - ^ - r * % On the first two mesh levels, 200 multigrid cycles are used. The fine mesh results 
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for a wing-body-winglet analysis at 25 and 50 multi-grid cycles are compared with 
the results at 200 cycles in Fig. 8. The differences in solutions are hardly 
discemable after 50 multigrid cycles on the fine mesh. To obtain a similarly 
converged solution without using successive mesh refinement requires more than 
i multigrid cycles in the fine mesh. For transonic cases with strong shocks, 

successive mesh refinement reduces the multigrid cycles from 400 to 200 on the 

fine mesh. The shock capturing scheme in the present Euler solver typically 
smears out the shock discontinuities over three to five grid points. The 
convergence of the shock discontinuities is more efficient when successive mesh 
refinement is used. 

For the consistency test with respect to dissipation parameters, Euler analyses 
for a wing-body-winglet configuration with two different sets of dissipation 
parameters were conducted. In Fig. 9, symbols show a solution using dissipation 
parameters twice as large as those for the solution given by the solid line. The two 
results are very close, although increasing the values of the dissipation parameters 
appears to smooth out some of the local variations in the solution. Important flow 
features, such as the peaks and locations of the leading edge expansions, are 
unchanged from the wing-root to the tip of the winglet. 

For the grid sensitivity test, the Euler solutions on two grids with significantly 
different grid quality were obtained for comparison. Both grids contain 200,000 
points and have identical surface grid distributions. However, Grid-1 has a strong 
i****>^m kink at the tip of the winglet, whereas Grid-2 is smooth near the tip of the winglet. 

The differences in solutions due to these grid effects are small and localized (Fig. 
10). Away from the tip of the winglet, the differences are essentially non- 
discemable. Again, important flow features such as shock locations are unchanged 
throughout the entire wing. 

” « A convergence study of the Euler solutions with respect to mesh refinement 

: ... - requires the use of an extremely fine mesh. This becomes prohibitively expensive 

in 3D. Therefore, a 2D version of the Euler code was employed for such studies. 
We chose the NLR7301 shockfree airfoil as a test case. Four successive levels of 
grids were used starting from grid-1 of 40x8, then grid-2 of 80x16, grid-3 of 
160x32, and ending with grid-4 of 320x64 (Fig. 11). The Euler solution in grid-1 
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is noticeably different from the finer grid solutions. The airfoil geometry is not 
very well defined in this grid. As the grid is refined, it is seen that the solutions on 
grid-2 and -3 are reasonably close to the finest grid solution on grid-4. This 
suggests that the 160x32 grid is adequate for this type of calculation. - : > •< 


4. Embedded Flow Solution 

A capability has been developed for the computation of the detailed flow 
characteristics in the local vicinity of an aft-mounted propfan engine in the 
presence of a flowfield induced by a complete aircraft configuration at subsonic 
and transonic speeds. The method used in the present study is to embed fine grids 
about the aft-mounted propfan engine and to solve the Euler equations on this grid 
using boundary conditions interpolated from a global Euler solution. Details on 
the local grid generation and the embedded procedure for the local flow field are 
discussed in the following subsections. 

a. Embedded-Region Grid Generation 

To resolve a detailed flow field in the propfan nacelle region, grids with suitable 
topology are generated. For a stmt-mounted, domed nacelle with propeller disk, 
shown in Fig. 12, a cylindrical grid that wraps around the nose of the nacelle is 
used. In the case of a propfan nacelle with inlet and exhaust, a cylindrical C-type 
grid wrapped around the fan cowl leading edge, shown in Fig. 13, is used. In 
either case, an H-type grid is used along the stmt. The external boundary is 
obtained by constmcting circular cross section in the streamwise direction with a 
given radius R, where R is set equal to the distance between nacelle centerline and 
the stmt/fuselage intersection. The upstream and downstream boundary can be set 
by the user through input. The number of grid points in each coordinate direction 
are also input by the user. A typical stmt-mounted, domed nacelle shown in Fig. 12 
uses 101x21x25 grids with 81x25 points on the nacelle surface, 81x21 points on 
the stmt, and 13x25 points on the propeller disk. The local grid is generated 
independently of the global configuration and the global grid system. The outer 
boundary surface of the local grid does not have to be one of the grid faces of the 
global system. Interpolations of the global flow variables onto the local grids are 



needed in order to establish the boundary condition as well as the initial condition 
for the local embedded analysis. 

b. Interpolation method 

The global flowfield solution can be computed as described in Section 3. A 
dense local grid about the aft-mounted propfan engine can then be generated. 
Such a grid will allow resolution of such geometry details as the inlet lip and strut 
leading edge. The solution from the global flowfield analysis are extracted and 
used as boundary conditions for the embedded region through interpolation (Figs. 
12, 14). More details about the interpolation method are given in Appendix D. 

c. Embedded Method 

An Euler solver based on the same technology used in the global solution 
method is developed to solve the detailed flow field around the aft-mounted 
propfan engine. The interpolated global flow solution is used as an initial 
condition for the embedded flowfield and the boundary condition in the far-field. 
As the solution evolves, the far-field boundary condition is updated according to 
characteristic relations derived from the local one-dimensional Euler equations 
shown in Appendix C. The implementation of a characteristic boundary condition 
is needed in order to prevent the disturbances reflecting back to the flowfield. 
Other types of boundary conditions are implemented to the solver in the same way 
as the global solution method and are described in Appendix C. 

The present grid embedding approach assumes that the global flow solution 
provides adequate resolution in the region away from the local embedded region, 
such that the interpolated data can be used as boundary condition for the embedded 
region. No iterative procedure between the global solution and the embedded 
solution is used in the present work. A better approach of resolving flowfield 
details around a complex aircraft configuration is to use a general multiblock 
approach, where a high resolution multiblock grid can be generated and used for 
each individual block. 
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5. Results and Discussions 


3 Results from global flowfield analysis will be presented first followed by the 
-*» embedded solution. How to use the program for these analyses is described in 
Volume II (Ref. 15). 

a. Global Flow Solution 

An aft-mounted propfan configuration shown in Fig. 15 is analyzed to 
demonstrate the capability of the out-of-core Euler flow solver. The grid system 
consists of 241x37x33 points with 81x25 points on the wing, 65x9 points on the 
stmt, 37x17 points on the horizontal tail, 27x13 points on the vertical tail, 24x69 
points on nacelle and plume, and 2925 points on the fuselage. 

The flow conditions are set at Moo = 0.8 , a = 1.5 deg. The propeller power 
conditions are thrust coefficient = 0.1232, zero normal force and side force 
coefficients. The isobars for the complete aircraft are shown in Fig. 16. Figures 
17 through 19 show the detailed Mach contours in the aftbody regions for 
propeller-on and propeller-off cases. It can be seen that the flow accelerates 
substantially in the aftbody regions for the propeller power-on case. Fig. 19 
reveals the formation of a supersonic zone on the lower surface near the nacelle- 
stmt junction. This occurs for both propeller power-on and power-off cases. The 
analysis identifies potential aerodynamic problems of the configuration such as 
wave drag penalty. 

The analysis information provides opportunities for engineers to recontour or 
redesign the problem areas of the aircraft. The present analysis capability can 
therefore complement the wind tunnel test to render a superior aerodynamic design 
in the aircraft development process. The effectiveness of the analysis capability 
developed in the present study can be further improved by coupling our capability 
with an automated iterative procedure for aerodynamic design developed at NASA 
Langley (Ref. 2). 

A crucial issue for aft-mounted propfan engine installations is the wing wake 
* location relative to propeller blade location. Specifically, the intersection of wing 
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wake with propeller blade should be avoided. The streamline tracing capability 
developed in the present study can serve this function very well. The streamlines 
i emitted from the trailing edge of the wing for the previously analyzed case are 
shown in Fig. 20. 

b. Embedded Flow Solution 

A stmt-mounted domed, nacelle shown in Fig. 12 is analyzed to demonstrate 
the embedded Euler solution. The far-field boundary condition is provided by the 
global flowfield solution through interpolation, and the propeller power conditions 
are the same as that in the global solution. Figures 21 and 22 show the Mach 
number contours of the embedded nacelle-stmt from the embedded Euler solution. 
This trend is similar to the results from the global solution; however, considerably 
more detail can be seen in the embedded solution (Fig. 23) due to the use of a high 
resolution local embedded grid. 

The grid embedding can provide a very useful capability for predicting the 
detail flow characteristics in the local vicinity of an aft-mounted propfan engine in 
the presence of a flowfield induced by a complete aircraft. However, the 
procedure for embedded flow solution (described in Section 4) involves several 
steps and would require that the global solution be carried out first. Furthermore, 
the current grid embedding procedure does not allow the local embedded solution 
to feed back to the global flowfield. To fully couple the global and local flowfield 
"would require additional iteration. A better approach of resolving flowfield details 
around a complex aircraft configuration is to use a general multiblock approach, 
where high resolution multiblock grid can be generated and used for each 
individual block. 


6. Conclusion 

A highly efficient and accurate Euler solver for predicting the 
airframe/propulsion integration effects for an aft-mounted turboprop transport 
aircraft configuration has been developed. The solution method is insensitive to 
grid quality and dissipation parameters. Convergence of the solution with respect 
to mesh refinement was demonstrated in 2D to further assess the numerical 
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accuracy of the Euler solutions. An embedded flow solution capability has been 
developed for predicting the detailed flow characteristics in the local vicinity of an 
aft-mounted propfan engine in the presence of a flowfield induced by a complete 
aircraft. 
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Appendix A 

Grid Generation Equations 


A grid generation method based on the solution of a set of elliptic grid 
generation equations 


^xx + ^yy + ^zz — P 

(Al) 

T|xx + %y + Bzz = Q 

(A2) 

Cxx + Cyy + Czz = R 

(A3) 

are used in the present work. Here x, y, z are physical coordinates and rj, £ are 
the computational coordinates. P, Q, and R are the added source terms to control 


the grid spacing in the flow field. The strength of the source terms is dictated by 
the boundary grid distributions and is written in terms of the first and second 
derivatives of the arc length distributions along the boundary surfaces by 


P = <Si)&/(SiX 

(A4) 

Q = '(Sj )r|ri / (Sj)r[ 

(A5) 

R = -(S k )cc / (Sk)c 

(A6) 


where Sj, Sj, S k represent arc length of a grid line along the boundary surface. 
These source terms are computed along each boundary surface and are interpolated 
into the flow field using a simple linear interpolation method. Because the source 
strength varies smoothly from one boundary surface to another, the interior grid 


distribution will also vary smoothly from one surface to another surface. 
Additional weighting functions can be introduced to further increase the 
controllability of the field grids. 
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Appendix B 

Numerical Method for the Euler Solver 

The basic numerical scheme follows Jameson’s finite volume formulation of 
solving the Euler equations 

dw/dt + df/dx + dg/dy + dh/dz = 0 (B 1) 


-p - 


- pu - 


- pv - 


- pw i 

pu 


P+pu 2 


pvu 


pwu 

pv 

f = 

puv 

g = 

P+pv 2 

h = 

pwv 

pw 


puw 


pvw 


P+pw 2 

-PE- 


- puH - 


- pvH - 


- pwH - 


Here, the conventional notations p, u, v, w, P, E are used for density, three 
components of velocity, pressure, and the total energy. The total enthalpy H is 
defined as 


H = E + P/p = Y P/(Y -l)p + (u 2 + v 2 + w 2 )/2 


(B2) 


with y representing the ratio of specific heat. In the governing equations, the 
density and the pressure are normalized with respect to its freestream values. The 
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total energy and the total enthalpy are normalized with respect to Poo/poo, and the 
velocity is normalized with respect to the square root of Poo/poc. 

The unsteady Euler equations are solved by the multistage Runge-Kutta time- 
stepping scheme 

w^ = w n 

w^ = w^ - cci AtR^ 

(B3) 

w (m-l) _ w (0) _ (X m _ 1 AtR ( ' m-2 ' ) 

W n+1 = w (m) = w (0) - AtR (m_1) 
with R (m) = C(w) (m) + D(w) (m) 


Here, C(w)^ m ^ represents convective flux terms, and D(w/ m - ) represents artificial 
dissipation terms added to control numerical oscillations. Central difference 
approximations are used for both convective and dissipative term calculations. A 
spectral radius scaling on the dissipation terms is used to minimize the amount of 
artificial dissipations, as shown in the following 


D(w)i = d i+ i/2,j,k - di-i/2,j,k (B4) 

d i+ i/2,j.k = W,j.k(e G) - 8 (4) 5^ 2 )(W i+ i i j ik - Wj,j. k ) (B5) 

where ^+\i 2 jx is the spectral radius at face i+1/2, and is defined as 

^i+l/2,j,k = (fln i+1 /2,j,k + a i+l/2,j,k)Si+l/2,j,k (P6) 
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Within each block of the flow field, a V-cycle multigrid scheme is used 
provided the number of cells in each coordinate direction can be divided by 2 to 
the nth power, with n+1 representing the number of multigrid levels. — 
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Appendix C 

Boundary Condition Implementations 


The boundary conditions encountered in the complete airplane analysis include: 
1) Solid-surface boundary condition, 2) Far-field boundary condition, 3) Nacelle- 
inlet boundary condition, and 4) nacelle exhaust plane/propeller disk boundary 
conditions. Each type of boundary condition implementation is described in the 
following subsections: 

1). Solid-surface boundary condition 

Along the configuration surface, the tangency condition is enforced by setting 
q • n = 0. For the cell-centered scheme used in the present work, the pressure on 
the configuration is also needed to evaluate the convective flux terms in the 
momentum and energy equations. Simple extrapolation from the flowfield solution 
can be used. However, to improve numerical accuracy, a normal momentum 
equation that relates the normal derivative of pressure to the streamwise and 
spanwise derivatives is used, i. e., along surface 1=2, one has 


s 2j Sij ap/a^i - pUiUja(s 2 j)/a^ = o 

with Sy = h a^i/axj , and Ui=SjjUj, h= lOxj/O^jl 


(Cl) 


2). Far-field boundary condition 

Because the calculation is performed in finite region, correct implementation of 
far-field boundary condition is essential for numerical accuracy as well as for 
numerical stability. In the present study, the characteristic boundary condition 
based on locally one-dimensional Euler equations is used. Along each face of the 
far-field boundary cell, the normal velocity is computed and checked. For 
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subsonic freestream condition, if it is an inflow boundary, one Riemann variable 
extracted from the interior flow field together with the tangential velocity, the 
enthalpy, the entropy and another Riemann variable evaluated from upstream 
values are used to determine all flow variables, i.e.. 


q t =q too , 8 = 800, H = Hoo , 

qn + 2a/(y - 1) = q^ + 2a e /(y - 1) , (C2) 

qn - 2a/(y - 1) = - 2aoo/(y- 1) . 


For the subsonic outflow boundary, one Riemann variable based on freestream 
values together with four extrapolated quantities are used to calculate the flow 
variables, i.e., 


q t = flte > S = S e , H = H e , 

qn + 2a/(y - 1) = q^ + 2aoo/(y - 1) , (C3) 

q„ - 2a/(y - 1) = q^ - 2a e /(y- 1) . 


with subscript e representing values extrapolated from the interior flowfield. In the 
case of supersonic inflow boundary, all variables are assigned to their freestream 
values. For supersonic outflow boundary, all flow variables are extrapolated from 
the interior flow solutions. 


3 ). Nacelle-inlet boundary conditions 

Along the nacelle inlet face, it is treated as an outflow boundary. Three choices 
of boundary conditions can be used: 
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a). Mass flux boundary condition 


If the capture stream tube area of the nacelle, i.e., A**, is given, the mass flux 
at the inlet face is given by 


(Pfln)f (Poo Aoo(joo)/Af 


(C4) 


with Af representing inlet face area, 

b). Velocity boundary condition 

For a given captured stream tube area A^ , one can convert the mass flux 
boundary condition into the velocity boundary condition using the isentropic 
relation to relate density to the local velocity, i.e.. 


qn/q„ = (A„/A f )[l - (Y-l)(M,») 2 {(q/q M ) 2 - (C5) 


A simple Newton’s iteration procedure is written to solve the above nonlinear 
equation for q„. 

c). Pressure boundary condition 

The static pressure P along the inlet face can be prescribed. The rest of the 
flow variables are extrapolated from the upstream flow field. 

4). Nacelle exhaust plane/propeller disk boundary conditions 

Along the nacelle exhaust plane, as well as along the downstream side of the 
propeller disk, similar boundary conditions can be applied. In either case, it is 
assumed to be a subsonic inflow boundary. The first type of boundary condition 
that can be applied to both nacelle and propeller disk prescribes the total pressure, 
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total temperature together with the swirl along the boundary. These quantities are 
then used to evaluate the flow variables according to the following formula 


M 2 = 5q 2 /(7(T t ) 2 -q 2 } , 

p = P t /(1 + 0.2M 2 ) 3 ' 5 , (C6) 

P = 7P/KY-l)(CpT t - q 2 /2} 


with q obtained from extrapolation. Another type of boundary condition 
frequently used for propeller simulation prescribes the thrust per unit area F x , the 
normal force per unit area F y , the side force per unit area F z , and the power done 
by the propeller Q to calculate the flow variables downstream of the propeller disk, 
i.e.. 


pUu + S X P - piUjUi + S x Pj + SF X , 
pUV + SyP = PlUlVj + SyP J + SFy , 

pUw + S Z P = PiUiWj + S z Pi + SF Z , (C7) 

u 2 + v 2 + w 2 = q 2 , 
p = 3.5P / {C p T tl + SQ/pU - q 2 /2) 


with subscript 1 denoting upstream conditions and S x (or S y , S z ) is the projection 
area in x- (or y-, z-) direction of a cell face on the propeller disk, and S is the cell 
face area. In either type of propeller disk boundary condition method, the 
boundary condition along the upstream side of the propeller disk is obtained by 
setting the mass flux to be equal to the mass flux of the downstream side. 
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Appendix D 

Interpolation Method 


• For the aft-mounted turboprop transport, the test case consists of 241x37x33 
grid points (240x36x32 cells), the flowfield is divided into 3x2 blocks. There are 
three blocks in the J (normal) direction and two blocks in the K (spanwise) 
direction. There is no blocking in the I (wrap-around) direction. The cell 
dimension of each block is 240x12x16. 

A reduced flowfield is obtained by using data from the inner two blocks in J 
direction. In addition, flow variables are calculated at the grid points through 
averaging. The reduced grid (and flowfield solution) consisting of 239x24x33 grid 
points is written in single block. The reduced flowfield will be used for 
interpolating the data onto the local embedded grid. 

The local embedded grid consists of 101x21x25 grid points (100x20x24 cells). 
The cell center locations are computed for all cells as interpolation points for the 
interpolation of flowfield data from the reduced global solution. A cell-centered 
local embedded grid is then written in single block. 

Within a cell in the reduced global flowfield, a flow variable can be interpolated 
using trilinear basis functions for any interpolation point within the cell, including 
the cell boundary. Let F be a variable (F may be a flow variable or one of the three 
physical coordinates), then within a cell 

8 

F= L F m N m (Dl) 

m=l 


where 


Fi = Fi,j, k 


24 



F 2 = Fi+l,j,k 
F 3 = F i,j+l,k 
F 4 = F i+l,j+l,k 
^5 = Fi,j,k+1 
f 6 = Fi+l,j,k+l 
F 7 = Fi,j+i,k+l 

f 8 = Fi+i.j+i.k+i 

and N m ’s are the trilinear local basis functions defined as: 

Ni -(1 -©(1-10(1-0 

N 2 =4(1 -T|)(l -C) 

n 3 = (1 - 4) ri (1 - 0 
N 4 = 4 r» (1 - 0 
n 5 =(1-4X1-11)4 
n 6 =4(1 -t i)C 
n 7 = (1 - 4) n C 
n 8 =4 t i C 

The 4 tj and C, are the local coordinates bounded from 0 to 1 
£ = 0 at i and % = 1 at i+1 
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r\ = 0 at j and r\ = 1 at j+1 

£ = 0 at k and £ = 1 at k+1 

Let a point P with physical coordinates (x, y, 
set F = x (or y, z) in Eq. (Dl) to obtain 


-x ■ 


■ x m ‘ 


8 1 


y 

II 

M 

Z 

3 

y m 


m=l 


- z - 


■ z m ‘ 


z) be in the cell for interpolation. We 


(D2) 


We then solve Eq. (D2) for the local coordinates £, r|, and £ through Newton 
iteration starting from cell center values, 

£ = 0.5 

r| =0.5 

£ = 0.5 

The computed local coordinates £, r|, and £ can then be used for interpolating a 
flow variable (e. g. pressure p) at the point P by setting F = p in Eq. (Dl). 

Given an interpolation point P(x, y, z) in the local embedded grid, we need to 
find the cell in the (reduced) global flowfield that contains this point. Let a cell in 
the global flowfield be referenced by its lowest i, j, k indexes. For example, the 
cell in the previous discussion will be referenced by i, j, k. We may begin the 
search procedure at the cell corresponding to i=l, j=l, and k=l. We can then use 
Eq. (D2) to compute the local coordinates £, T) and £. When any of the local 
coordinates is out of bound, the point P is considered not in the current cell. The 
values of the local coordinates will be used for stencil walk. 
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If ^ < 0 then move to a cell at i-1 
If £ > 1 then move to a cell at i+1 
If r\ < 0 then move to a cell at j-1 
If r| > 1 then move to a cell at j+1 
If £ < 0 then move to a cell at k-1 
If C, > 1 then move to a cell at k+1 

Increment (or decrement) of the indexes will not take place if the resulting cell is 
not within the reduced global flowfield. 
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Fig. 1. An Advanced Propfan Airplane Configuration 
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Fig. 2. Grid Generation and Blocking Strategy 
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Fig. 3. Grid Topology at a Given Span Station 
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Fig. 4. Field Grids on Spanwise Cuts 
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Fig. 5. Field Grids on Streamwise Cuts 
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WING/BODY INTERSECTION AFTBODY SECTION 


Fig. 6. Grids on the Wake 
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Fig. 7. Surface Grids and Partial Field Grid for an 
Aft-Mounted Propfan Airplane 
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Fig .8. Convergence in Multigrid Cycles 
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Fig. 13. C-Grid for Nacelle with Inlet and Exhaust 
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Fig. 14. Complete Flo 
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Fig. 15. Surafce Grid for 
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Isobars for an Advanced Propfan Aircraft at Macb=0.8, Alpha=1.5 
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acli Contours on Fuselage Aftbody/ Vertical Tail for an Advanced 
opfan Configuration at Mach =0.8, Alpha =1.5 
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Fig. 18. Mach Contours in Aftbody Regions for an Advanced Propfan 
Configuration at Mach =0.8, Alpha =1.5, Upper Surface 





19. Mach Contours in Aftbody Regions for an Advanced Prop fan 
Configuration at Mach =0.8, Alpha =1.5, Lower Surface 













Fig. 20. Streamline Tracing for an Aft-Mounted Propfan Aircraft 
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Comparisons of Embedded Solutions in the Propfan Engine 
Regions, Mach =0.8, Alpha =1.5, Upper Surface 


Global Solution 
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Fig* 2S. Comparisons of Global/Embedded Solutions in the Propfan 
Engine Regions at Mach =0.8, Alpha=1.5, Lower Surface 
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